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Abstract — We present a derivation of the flavor neutrino states which describe neutrinos produced 
or detected in charged-current weak-interaction processes, including those operating in neutrino oscil- 
lation experiments. We also present a covariant derivation of the probability of neutrino oscillations 
which is consistent with the fact that flavor is Lorentz-invariant. Finally, we clarify the negative an- 
swers to three commonly asked questions: "Do charged leptons oscillate?"; "Is the standard phase 
wrong by a factor of 2?" "Are flavor neutrinos described by Fock states?". 

1 Introduction 

The study of neutrino oscillations is a very important field of contemporary experimental ant theoretical research 
in high-energy physics. The main reason is that neutrino oscillations is a consequence of the existence of neutrino 
masses, as was discovered by Bruno Pontecorvo in the late 50's [l|, |2j, |3[ . However, neutrinos are massless in the 
Glashow-Salam- Weinberg Standard Model Hence, neutrino oscillations represents an open window on 

the physics beyond the Standard Model. 

In view of such importance of neutrino oscillations, it is necessary to have a correct and clear understanding of all 
aspects of the theory of neutrino oscillations, in order to have a correct interpretation of ex perimental results. The 
theory of neutrino oscillations has been discussed in ma ny pape rs I^l^l^ llfl[lllj l2 i ni 
ESQ EE 00 III 013 1110 and reviewed injplli 

(conference contributions and papers with other points of view are listed in Ref. |45]). 

In this report we first review briefly the standard plane-wave theory of neutrino oscillations in Section |2] In 
Section[3]we discuss the problem of the definition of flavor neutrino states and in Section|4]we show that the flavor 
neutrino states are appropriate for the description of neutrinos produced or detected in charged-current weak- 
interaction processes. In Section[5]we present a covariant derivation of the probability of neutrino oscillations. In 
Section[6]we try to give clear answers to three questions which are often asked. Finally, we conclude in Sectional 




2 Standard Plane- Wave Theory 

Neutrino oscillations are a consequence of neutrino mixing: the left-handed flavor neutrino fields v a L, with a = 
e, /I, t, are unitary linear combinations of the massive neutrino fields 1 ) VkL, with fc = 1,2, 3, 

3 

VaL = ^2 U ak VkL (« = e, fl, t) , (2.1) 

k=l 

where U is the unitary mixing matrix. 

In the standard plane-wave theory of neutrino oscillations 0,0,0, 03 E3l it is assumed that the neutrinos created 

1 ^For si mplicity we consider the simplest case of three-neutrino mixing, neglecting the possible existence of additional sterile neutrino fields 
(see Refs. SElHEI!!)- 



or detected together with a charged antilepton 2 ) £+, with a = e,fi, t, are described by the flavor neutrino state 
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Since the massive neutrino states \v k ) have definite mass m k and definite energy E k , they evolve in time as plane 
waves: 

i^Mt))=J?\Mt))=EkWk(t)) => \v k {t))=e- iE * t \v k ), (2.3) 
at 

where ffi is the Hamiltonian operator and \v k ) = \v k (t = 0)). The consequent time evolution of the flavor 
neutrino state ( 12. 21 is given by 

KW) = E^ e " <£ki W= E (E^e-^^J m, (2.4) 

fe=l 0=e,n,T \fc=l / 

which shows that if the mixing matrix U is different from unity (i.e. if there is neutrino mixing), for t > the state 
Wa(t)), which has pure flavor a at t — 0, is a superposition of different flavors. The quantity in round parentheses 
in Eq. ( 12. 41 is the amplitude of v a — ► z/g transitions at the time i after v a production, whose squared absolute value 
gives the probability of a v a — » z/g transitions: 



p„ a Wf) = K^K(*)>l 2 



fe=i 



2 3 



£ E£ fc t/^ C/ Qj e-'C**-^)* . (2.5) 



One can see that P 1Ja ^ L , fj (t) depends on the energy differences E k — Ej. In the standard theory of neutrino 
oscillations it is assumed that all massive neutrinos have the same momentum p. Since detectable neutrinos are 
ultrarelativistic 3 ), we have 

2# J 2£ 



Am 

£ fc = X /P z + ml ~ £ + ^ E k -E J = , (2.6) 



where E = |p| is the energy of a massless neutrino (or, in other words, the neutrino energy in the massless approx- 
imation), and Am|j = m| — m 2 . In order to obtain a measurable flavor transition probability, it is necessary to 
convert the time t, which is not measured in neutrino oscillation experiments, in the known source-detector dis- 
tance L. Considering ultrarelativistic neutrinos, we have t ~ L, leading to the standard formula for the oscillation 
probability 

^ / A 9 r\ 

P v ^ Vf) (L) = J2 Kk Up k U aj U* 0J exp -i _^|_ , (2.7) 

which is used in the analyses of the experimental data on neutrino oscillations in vacuum. 
Summarizing, there are three main assumptions in the standard theory of neutrino oscillations: 

(Al) Neutrinos produced or detected in charged-current weak interaction processes are described by the flavor 
states i2.2\. 



(A2) The massive neutrino states \v k ) in Eq. ( 12. 21 have the same momentum. 

(A3) The propagation time is equal to the distance L traveled by the neutrino between production and detection. 



3 Flavor Neutrino States 

The flavor neutrino states which describe neutrinos produced or detected in charged-current weak interaction pro- 
cesses can be derived in the framework of Quantum Field Theory i32ll36ll42Tl . The final state \ f) in an interaction 
process is given by 

\f) = S\i), (3.1) 

2 ) Since a is an index, we use the notation i e = e, £ M = fi, l T = r. 

3 ) lt is known that neutrino masses are smaller than about one eV (see Refs. l48tl5Cll V Since only neutrinos with energy larger than about 100 
keV can be detected (see the discussion in Ref. I32l0 . in oscillation experiments neutrinos are always ultrarelativistic. 
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where \i) is the initial state and S is the S-matrix operator. In general, the final state can be a superposition of 
orthogonal states, 



\.f) = Ak ' ( 3 - 2 ) 

with (fk\fj) = 5kj- From Eq. (13 . 1 1 it follows that the amplitudes Ak are given by 

-4* = </ fc |/> = </ fc | fif |*> . (3.3) 

Turning now to the case of neutrinos, for definiteness let us consider neutrino production in the general decay 
process 

Pi ^ P F +e+ + is a , (3.4) 

where Pj is the decaying particle, Pp is a decay product (which can be absent, for example in tt + — > /i + + z/„), 
and a — e, /i, r. Other production processes, as well as detection processes, can be treated in a similar way with 
straightforward modifications to the formalism. The final state of Pj decay, 

\f)=S\Pj), (3.5) 

is given by 

\f) = N a \v a ,t+,P F ) + ... , (3.6) 

where N a is a normalization coefficient and the dots represents other decay channels and the undecayed initial 
state Pj itself. The state \v at £+, Pp) can be decomposed in a superposition of orthogonal states corresponding to 
different massive neutrinos: 

1 3 

K,l+,P F ) = j r Y,A ak \v k ,et,P F ), (3.7) 
a fe=i 

with the amplitudes A a k given by 

A ak = {u k ,l+,P F \f) = (v k J+,P F \S\Pi). (3.8) 
The normalization of the state (13.71 requires that 



\N a \ 2 = ^\A ak \ 2 . (3.9) 



fe=i 



It is convenient to write Eq. d3.7t in the form 



/ 3 \ _1 / 2 3 

\P F Jt)K) = \P F ,£i) El- 4 ^! 2 J2-A ak \v k ), (3.10) 



\fe=l / fe=l 

which shows that the final neutrino in the decay d3.4i is described by the flavor state 

/ 3 \ _1 / 2 3 

K)=[J2\A ak \ 2 \ ^^ Qfc |^ fe ), (3.11) 

which is calculable through Eq. J3.8I >. This flavor state has the same structure as the standard one in Eq. d2.2l >. 

/ 3 \-Va 

In order to see if the amplitudes ^4 a fc I |«4 Q fe | 2 ] agree with the standard ones given by U* k , we need to 

\fc= i / 

calculate the matrix elements in Eq. d3.8l >. In the case of the low-energy weak decay ( 13. 41 the effective S-matrix 
operator can be approximated at first order in the perturbative expansion as 

S ~ 1 - i J d i x,^ I cc (x) , (3.12) 

where J^_ cc (x) is the charged-current weak interaction Hamiltonian 

M{ cc (x) = ^ J2 ^)/(l-7 S )4W^ ft W+H,c. 

" a—e,fi,r 



3 



G 



^ £ E c/ ^^)^( 1 -^ 5 )^( a; ) J p p ^ PF ( a; )+ H - c - 



(3.13) 



a=e,/j,,r k=l 

The operator Jp^ PF (x) describes the transition Pj — ► Pp. For the matrix element (uk, l^i Pf\ S \Pj) we obtain 



with 



Gy 
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Aak = (v k ,tt,PF\S\P I )=U* k M ak , 



d 4 x (v k , C, Pf\T^(x) Y (1 - 7 5 ) i a {x) J^ Pp (x) |Pr) 



Hence, the flavor neutrino state J3 .111 is given by 



■'..> = £■ 



Mak 



k=l 



Kk Wk) 



(3.14) 
(3.15) 

(3.16) 



We notice that the flavor neutrino state J3.16> has the structure of the standard flavor state 12.21 . with the amplitudes 
of the massive neutrinos \ vk) proportional to the corresponding element U* k of the mixing matrix. The additional 

-1/2 



disappear in the realistic case of ultrarelativistic neutrinos, because 



coefficients M a k £ \ U a j\ 2 \M a j 

in this case the experiments are not sensitive to the dependence of M a k on the neutrino masses m k , leading to the 
approximation 

Mak -M a => — =-1, (3.17) 



^2\U ao f\M c 

\ 3 = 1 



where we have used the unitarity relation | t/a j | 2 = 1. 



3=1 

Hence, in the case of ultrarelativistic neutrinos we obtain the standard flavor neutrino states (12.21 in a rigorous 
way in the framework of Quantum Field Theory. Since in all neutrino oscillation experiments the approximation 
d3 . 171 is valid, in neutrino oscillations flavor neutrinos are correctly described by the standard flavor states ( 12. 21 , 
confirming the validity of the standard assumption |(Al)| Let us emphasize that the validity of this approximation 
is very important, because the standard flavor states (12.21 do not depend on the kinematics of the process in which 
the flavor neutrino is created or detected. Hence, it is possible to derive a general expression for the oscillation 
probability, as we will see in Section|5] 



4 Neutrino Production and Detection 

In this Section we show that the flavor states ( 13. 1H are appropriate for the description of neutrinos produced or 
detected in charged-current weak interaction processes. We consider for simplicity only the decay process (13. 4L 
but the same reasoning can be applied with straightforward modifications to any production or detection process. 

It is well known that the probability of the decay (13. 4> is the sum of the decay probabilities 

in the different massive neutrinos v k weighted by the squared absolute value |J7 a fe| 2 of the element of the mixing 
matrix that weights the contribution of v k to the charged-current weak interaction Hamiltonian ( I3.13l l. The reason 
is that the massive neutrinos have definite kinematical properties and constitute the possible orthogonal asymptotic 
states of the decay. In other words, each decay in a massive neutrino constitutes a possible decay channel. Hence, 
the probability of the decay ( I3.4t is given by 

3 

\A a ? = Y J \ lJ o,k\ 2 \MakW (4.1) 

k=l 

with Mak given by Eq. ( I3.15l l. 
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Let us now check that the description of the neutrino v a produced in the decay J3.4I > through the flavor state ( 13. Ill 
gives the correct result (14. II . Indeed, the amplitude of the decay process P.4t is given by 

/ 3 \ - 1 / 2 3 / 3 \ 1 / 2 

A a = (v a ,i>+,P F \S\P I )= £|A*fe| 2 J2 A -ki^,li,PF\S\P I )= ^|^ Qfc | 2 ■ (4.2) 

\fc=l / fc=l \/c=l / 

Therefore, the decay probability is given by an incoherent sum of the probabilities of production of different 
massive neutrinos: using Eq. (13. 14i we obtain 

3 

\A a \ 2 = J2 \^ ak \ 2 = \U ak \ 2 \M ak \ 2 , (4.3) 
fe=i 

in agreement with Eq. ( 14. 1> . Thus, the coherent character of the flavor state ( 13. Hi is correctly irrelevant for the 
decay rate. 

It is important to understand that this result is very important for the interpretation of the data of neutrino oscillation 
experiments, because in the analysis of these data it is necessary not only to calculate the oscillation probability 
but also the neutrino production and detection rates. We have shown that both tasks can be accomplished in a 
consistent framework with the description of the produced and detected neutrinos through the flavor states ( 13.111 1. 



5 Covariant Plane Wave Neutrino Oscillations 

Since flavor is a Lorentz-invariant quantity (for example, an electron is an electron in any reference frame), the 
probability of neutrino oscillations is Lorentz invariant |29, 35]. Indeed, the standard oscillation probability for- 
mula i2.7\ is Lorentz-invariant, as shown in Ref. [35]. However, the standard derivation of the oscillation proba- 
bility reviewed in Section|2]is not formulated in a covariant way. In this Section we present a covariant derivation 
of the oscillation probability which shows explicitly its Lorentz-invariance ll29tl3ttl42u56tl57l.l5o11 . 

In the following derivation we consider ultrarelativistic neutrinos. Of course, neutrinos are not ultrarelativistic in 
all reference frames 4 ) . However, since the oscillation probability is Lorentz-invariant, it can be evaluated in any 
frame and we can choose to evaluate it in a frame where the massive neutrinos are ultrarelativistic. Let us emphasize 
that this is not a peculiar choice, because there is a continuous infinite set of reference frames in which massive 
neutrinos are ultrarelativistic, which includes the one in which the detector is at rest. Indeed, we will always 
work with explicitly Lorentz-invariant expressions for the oscillation probability, because these expressions must 
be valid in the continuous infinite set of reference frames in which massive neutrinos are ultrarelativistic. 

The amplitude of v a — > vp transitions at a distance L and after a time T from the production of the flavor neutrino 
v a is given by 

A aP {L,T) = ( Vf} \e- lET+lPL \ Va ), (5.1) 

where E and P are, respectively, the energy and momentum operators. Since the massive neutrinos have definite 
masses and kinematical properties, using the flavor states i2.2\ , which are valid for ultrarelativistic neutrinos, we 
obtain 

3 

A a0 (L, T) = £ K k U [jk e~^ T +^ L , (5.2) 
fe=i 

with 

E k = ^pt + ml. (5.3) 

The oscillation amplitude ( 15. 2i is explicitly Lorentz-invariant. Let us notice that we do not adopt the standard 
assumption |(A2)| of equal momentum for different massive neutrinos. Indeed, it can be shown that in general 
massive neutrinos do not have neither equal momenta nor equal energies Il29ll30ll33ll33ll . However, in spite of the 
unrealistic character of the standard assumption |( A2)| we will see that the standard expression for the oscillation 
probability is correct. 

In oscillation experiments the neutrino propagation time T is not measured. As in the standard plane-wave theory 
reviewed in Section[2] for ultrarelativistic neutrinos it is possible to approximate T ~ L (assumption| (A3)t in the 
phase in Eq. \5.2i . The reason is that in reality neutrinos are described by wave packets Jlllfl2lfl4lll(Tlll8l[l9lE3l 

4 'We are grateful to Prof. D. V. Ahluwalia-Khalilova for asking by e-mail a stimulating question on this point. 
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us , which are localized on the production process at the production time and propagate 
between the production and detection processes with a group velocity which is close to the velocity of light. If 
the massive neutrinos are ultrarelativistic and contribute coherently to the detection process, their wave packets 
overlap with the detection process for an interval of time [T — AT , T + AT], with 

T = = * L ( 1 + V AT ~ a x , (5.4) 



v 



2E 2 



where v is the average group velocity, m 2 is the average of the squared neutrino masses, and a x is given by the 
spatial uncertainties of the production and detection processes summed in quadrature |25J (the spatial uncertainty 
of the production process determines the size of the massive neutrino wave packets). The correction Lm 2 /2E 2 
to T = L in Eq. 15.41 can be neglected, because it gives corrections to the oscillation phases which are of higher 
order in the very small ratios m 2 /E 2 . The corrections due to AT ~ a x are also negligible, because in all realistic 
experiments a x is much smaller than the oscillation length = AnE / Am\j, otherwise oscillations could not 
be observed I12l Il4 l33l 14311 . One can summarize these arguments by saying that the approximation T ~ L is 
correct because the phase of the oscillations is practically constant over the interval of time in which the massive 
neutrino wave packets overlap with the detection process. 



Using the approximation T ~ L the phase in Eq. (15.21 becomes 

where E is the neutrino energy neglecting mass contributions. Let us emphasize that both the left and right hand 
sides of Eq. (15.51 are Lorentz-invariant |35|. 



Inserting the expression (I5.5i in Eq. ( 15. 21 one can obtain in a straightforward way a probability of v a — ► vp 
transitions in space which coincides with the standard one in Eq. ( 12.71 . Therefore, the result of the covariant 
derivation of the oscillation probability confirms the correctness of the standard formula in Eq. ( 12.7b . 



It is important to notice that Eq. (I5.5i shows that the phases of massive neutrinos relevant for the oscillations are 
independent from the particular values of the energies and momenta of different massive neutrinos 0000 
15^1 . because of the relativistic dispersion relation (15.31 . Hence, as remarked above, the standard assumption |(A2)| of 
equal momenta for different massive neutrinos, albeit unrealistic, is irrelevant for the correctness of the oscillation 
probability. 



6 Questions 

In this section we try to clarify the negative answers to three questions which are often asked. 



6.1 Do Charged Leptons Oscillate? No! 

In order to understand the negative answer to this question it is important to comprehend that the only characteristic 
which distinguishes different charged leptons is their mass. Thus, flavor is distinguished through the mass of the 
corresponding charged lepton or, in other words, flavor and mass coincide for charged leptons. Since for charged 
leptons there is no mismatch between flavor and mass, there cannot be flavor oscillations of charged leptons 1 29] . 

InRefs. 00EH 

it has been claimed that the probability to detect a charged lepton oscillates in space-time. This 
claim has been refuted in Ref . 1 22] . A similar effect, called "Lambda oscillations", which has been claimed to exist 
l62H63ll for the A's produced together with a neutral kaon, as in the process 7r~ + p — ► A + K°, has been refuted 
in Refs. 1 261 16411 . Considering the pion decay process 7r~ — * fi~ + v^, the authors of Refs. 1 59, 60, 61, 65] argued 
that, since the final state of the muon and antineutrino is entangled, if the probability to detect the antineutrino 
oscillates as a function of distance, also the probability to detect the muon must oscillate. However, it is well 
known that the probability to detect the antineutrino, irrespective of its flavor, does not oscillate. This property is 
usually called "conservation of probability" or "unitarity" and is represented mathematically by the general relation 
~^2g Pv^yp — 1. Hence, this argument actually proves that the probability to detect a charged lepton does not 
oscillates in space-time! 
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6.2 Is the Standard Phase Wrong by a Factor of 2? No! 

It as been claimed that the standard phase of neutrino oscillations in Eq. 12.71 . 

Ami L 

*w = -hJ-» (fU) 

is wrong by a factor of two 1 66] or there is an ambiguity by a factor of two in the oscillation phase i67ll68l | 69l.l70l| . A 
similar discrepancy by a factor of two in the phase of K° — K° oscillations has been claimed in Refs. 16211711172 
and an ambi guity by a factor of two has been claimed in Ref. 1 73 ] . These claims, which have been refuted in 
Refs. iHlMHl64ll7l.l75il7^l7ll. stem from the foUowing fallacious reasoning. 

Different massive neutrinos propagate with different velocities 

= -=r , (6-2) 

where Ek and pk are, respectively, the energy and momentum of the neutrino with mass mk, related by the 
relativistic dispersion relation ( 15. 3> . According to the fallacious reasoning, the phases of the different massive 
neutrinos wave functions after a propagation distance L should take into account the different times of propagation 
of different massive neutrinos: 

i> fc = Pk L- E k t k . (6.3) 

The propagation times are given by 

L E k 

t k = — = —L, (6.4) 

Vk Pk 

which lead, in the relativistic approximation, to the phase difference 

Am? L 

A$ kj = 4* - *j f— ■ (6.5) 



This phase difference is twice of the standard one in Eq. ( 16.1 



Let us notice that in Eq. d6.3i we have considered the possibility of different energies and momenta for different 
massive neutrino wave functions, as we have done in the covariant derivation of the neutrino oscillation probability 
in Section|5] The authors of Refs. claimed that a correct way to obtain the standard oscillation phase 

is to assume the same ener gy fo r the different massive neutrino wave functions. Apart from the fact that this is an 
unphysical assumption 1 29, 30L l33ll33ll . it is not true that the disagreement of a factor of two disappears assuming 
the same energy for the different massive neutrino wave functions, as clearly shown by the above derivation of 
Eq. d6.5J >. in which the energies of the different massive neutrino wave functions could have been taken to be equal 
Indeed, even if the different massive neutrino wave functions have the same energy, the time contribution 
—Etk + Etj to the phase difference A<tkj does not disappear if tk ^ tj. This contribution has been missed in 

Refs. ISFzaFzl. 



The mistake in the fallacious reasoning described above is due to a wrong use of the group velocity (16. 2> in the 
phase ( 16. 31 . which cannot depend on the group velocity. The group velocity ( 16. 21 is the velocity of the factor 
which modulates the amplitude of the wave packet of the corresponding massive neutrino. In neutrino oscillation 
experiments the envelopes of the wave packets of the different massive neutrinos propagate with the corresponding 
group velocity and take different times to cover the distance between the source and the detector. However, these 
different propagation times have no effect on the phases of the wave functions of the different massive neutrinos, 
whose interference generates the oscillations. Since the interference is a local effectjk must be calculated at the 
same time, as well as in the same point, for all the massive neutrino contributions 129113 lll34l l64l 174 l75l 176117711 . as 
we have done in Section|5] Since only the amplitude of the wave function of each massive neutrino is determined 
by the wave packet envelope, the only possible effect of the different arrival times of the envelopes of the wave 
packets of different massive neutrinos is a reduction of the overla p of the different wave packets, leading to a 



it massive neutrinos is a reduction or tne overlap or tne 
decoherence effect . For more details, see Ref. 13 111 



6.3 Are Flavor Neutrinos Described by Fock States? No! 

It must be noted that the flavor state J2.2I > is not a quantum of the flavor field v a fl5ll . Indeed, one can easily check 
that the flavor state J2.2i is not annihilated by the flavor field v a if the neutrino masses are taken into account. It 
is, however, possible to construct a Fock space for the quantized flavor fields 

Jz3 EJ El il E3 Hi 

Then, 
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it is natural to ask if the flavor Fock states describe real neutrinos, produced and detected in charged-current weak 
interaction processes. Let us notice that, since the flavor Fock states are different from the standard flavor states 
&2.2\ and the "exact" flavor states (13 . 1 1 i that we have derived in the framework of Quantum Field Theory, if the 
answer to the question under discussion were positive, the theory of neutrino oscillations would have to be revised, 
as claimed in Refs. f79l lifll Isil li^l 1831 l84l l85ll . However, the answer is negative, as explained in Ref. to 
which we direct the interested reader for the detailed proof. Here we mention only that the unphysical character 
of the flavor Fock states can be proved by reductio ad absurdum: the description of neutrinos created or detected 
in charged-current weak interaction processes through the flavor Fock states leads to unphysical results. Let us 
emphasize that this fact precludes the description of neutrinos in oscillation experiments through the flavor Fock 
states, because these neutrinos must be produced and detected in weak interaction processes. Hence, the flavor 
Fock states are ingenious mathematical constructs without relevance for the description of real neutrinos. 

7 Conclusions 

We have shown that the standard flavor states \2.2\ used in the derivation of neutrino oscillations can be derived in 
the framework of Quantum Field Theory in the realistic case of ultrarelativistic neutrinos (Section[3j. In Section^ 
we have shown that the "exact" flavor states (13 . 1 1 i are appropriate for the description of neutrinos produced or 
detected in charged-current weak-interaction processes, taking into account the neutrino masses. In Section[5]we 
have presented a covariant derivation of the probability of neutrino oscillations which is consistent with the fact that 
flavor is Lorentz-invariant. Finally, in Section[6]we have discussed the negative answers to three questions which 
are often asked. In conclusion, we would like to emphasize that the standard expression \2.1\ for the probability of 
neutrino oscillations in vacuum has been proved to be correct in the framework of Quantum Field Theory and all 
contrary claims stem from some misunderstanding. 
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